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Abstract 

We have studied charmless two-body B decays involving i] and ?/ in the framework 
of QCD improved factorization appraoch. The spectator hard scattering mechanism for 
ij' production have been re-examined and extended, which contributions are incorporated 
consistently into the factorizable leading contributions. It is found that the conventional 
mechanism would give Br{B — > rj' K) 30 x 10^^ which agrees with predictions based 
on naive factorization approaches. With the incorporation of spectator hard scatter- 
ing mechanism contributions, theoretical predictions are improved much and the bulk of 
Br[B rj'K) are accommodated in the reasonable parameter space. We have also pre- 
sented calculations of g*g* — i]^'^ transition form factors within the standard hard scattering 
approach. It is shown that the new contributions in the modes such as B ^ rj' + vector 
and B — > rj + ir, p, K^*^ are small. Direct CP- violation in those decay modes are predicted. 
It is shown that the prospects for measuring direct CP-violation effects in B^ —>■ r/'K^, 
r/'vr^, rjir^, and 7]K^ are promising at B factories in the near future. 

PACS Numbers: 13.25Hw,12.15Hh, 12.38Bx 
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1 Introduction 



The first evidence of strong penguin was observed by CLEO[|T| in 1997 with the announcement 
of 

Br{B- 7]'K~) = {Q.5ti,l ± 0.9) x 10~^ 
fir(fi° ^ r]'K'^) = {A.7+1-1 ± 0.9) x 10-^ (1) 

which are 2 ~ 4 times larger than any theoretical predictions existed at that time. Driven by 
the unexpected large data, these decays modes have been investigated extensively 0, ^ ^ ||. 
As a result of the investigations in past years, the contribution of the conventional mechanism 
estimated by using the naive factorization and effective Hamiltonian for B decays would account 
for l/4~ 1/2 of the data. Some new mechanisms are proposed to explain the unexpected large 
rates of B decays to Kt]'. Namely 

A) large intrinsic charm content of 77' through the chain b ccs rj's or through b — >■ 
CCS sg*g* rj's^; 

B) strong penguin 6 — > enhanced by new physics; 

C) spectator hard-scattering mechanism through g*g* 77' |^. 

For type A mechanism, its magnitude is characterized by the parameter f^, defined by the 
matrix < ?7'|c7^75c|0 >= —if^ip^- To account for the large branching ratio of B ^ v'^i 
would be as large as (50 ~ 18O)Mel^|0]. However, the later analyses have shown as small 
as few MeV [0, P|. It is also realized that strength of the process b — > ccs — > sg*g* rj's is 
very small |T0[| . Generally compared with uncertainties in form factors and light quark masses in 
the estimations of nonleptonic B decays, the contribution of type A mechanism for 77' exclusive 
production is probably safe to be neglected. For type B mechanism, it would be very interesting 
to find signals of new physics beyond the standard model(SM) in these decays if the SM is indeed 
incapable of accommodating experimental data. In this paper, we study those processes in the 
SM and investigate the possibility whether the SM can accommodate the present experimental 
data with new approach for the hadronic dynamics of B decays. 

For type C mechanism, it may be promising. Unfortunately it depends on some unknown 
quantities: the transition form factor of Tj' —g*g* and the light cone distribution amplitudes(DA) 
of the mesons in the process. The prediction in Ref.0 should be improved and incorporated 
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with the predictions of the basic mechanisms consistently. 

While the data of these decays reported in 1997[|l| is still puzzlingly large for theorists, 
robust experimental investigations are in progress. Recently, using the full CLEO II/II.V data, 
CLEO Collaborations |ll|] have improved their previous measurements oi B ^ rj'K 



B{B+ r]'K+) = (80^^° ± 7) x 10-^ 

^ = {89tll ± 9) X 10-^ (2) 



with observations of two new decay modes 



B{B+ r]K*+) = (26.4+1^ ± 3.3) x 10-^ 

^ r/Jr°) = (13^1;^ ± 1.6) x 10"^ (3) 



and upper limits for other 12 decay modes involving rj or rj'. 



Theoretically, the importance and generality of the pioneer works of Polizer and Wise|12 



and factorization arguments of Bjorken|jIB| are renewed by Beneke, Buchalla, Neubert and 



Sachrajda with the formation of "QCD factorization" for B hadronic decays|Tj, |T^. The 
factorization formula incorporates elements of the naive factorization approach (as the leading 
contribution) and the hard-scattering approach( as subleading corrections), which allows us to 
calculate systematically radiative (subleading nonfactorizable) corrections to naive factorization 
for B exclusive nonleptonic decays. An important product of the formula is that the strong 
final-state interaction phases are calculable from the first principle which arise from the hard- 
scattering kernel and hence process dependent. The strong phases are very important for 



studying CP violation in B decays. Detailed proofs and arguments could be found in [15]. Here 
we recall briefly the essence of the QCD factorization formula as follows. 

The amplitude of B decays to two light mesons, say Mi and M2, is obtained through the 
hadronic matrix element {Mi{pi)M2{p2)\Oi\B{p)) , here Mi denotes the final meson that picks 
up the light spectator quark in the B meson, and M2 is the other meson which is composed 
of the quarks produced from the weak decay point of b quark. Since the quark pair, forming 
M2, is ejected from the decay point of b quark carrying the large energy of order of m^, soft 
gluons with the momentum of order of Aqcd decouple from it at leading order of AgcD/^b in 
the heavy quark limit. As a consequence any interaction between the quarks of M2 and the 
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quarks out of M2 is hard at leading power in the heavy quark expansion. On the other hand, 
the hght spectator quark carries the momentum of the order of Aqcd, and is softly transferred 
into Ml unless it undergoes a hard interaction. Any soft interaction between the spectator 
quark and other constituents in B and Mi can be absorbed into the transition formfactor of 
B —>■ Ml which could be extracted from semileptonic decays B Milu. The non-factorizable 
contribution to B ^ A/1M2 can be calculated through the diagrams in Fig.l, which turns out 
to be subleading oder corrections to factorizable amplitudes. 

In this paper we study B — > r]^'^M(M = 'K,K^*\p) decays within the framework of QCD 
factorization approach |]T^, |15[. We compare our numerical results with the experimental data 
presented by CLEO collaboration []TT| . We find that the conventional mechanism contributions 
to B^ — i> ri'K~^ and B^ rj'K^ are about 27 x 10^^ and 28 x 10"^ respectively, which agree 
with many theoretical expectations based on naive factorizations. To explain the experimental 
data, contributions from new mechanisms with the strength as large as 40% ~ 50% of the 
strength of the conventional mechanisms are found. Incorporating the contribution of spec- 
tator hard-scattering mechanism(SHSM) to these decays, the experimental data could be well 
accommodated in the SM. SHSM is found to be important for B — * rj'K but not for B rj'K*. 

Our predictions agree with the data of measured branching ratios or lie below their upper 
limits of other decay modes. We also give our predictions of direct CP asymmetries in these 
decay modes. Direct CP violations in the four observed B decay modes B^ K~^*r] and 
B^ K^ri'[K^*ri) are found to be about a few percentages, but it can reach 10% for B^ — >• 
K^rj' . Considering its large branching ratio, we may expect that direct CP violation effects 
in charmless B decays would be firstly observed in B^ K^r]'. Large direct CP violation 
asymmetries are predicted for the decay modes B"^ v''^^^ and B^ r]K~^. Prospects 
for observing direct CP-violations in the these decay modes at B factories are very promising. 

This paper is organized as following. In Sec. 2 we present notations and calculations of the 
conventional mechanisms contributions to these decays. In Sec. 3 we calculate g*g*—ri' transition 



form factors using the standard hard scattering framework of Brodsky and Lepage pq]. With 
the form factor, we re-examine the contribution of the spectator hard scattering mechanism for 
B — s> ?7'-ft'0 and generalize it to other 14 decay modes. Section 4 contains our numerical results 
for the branching ratios of two body charmless B decays involving rj and t]'. Direct CP- violations 
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in these decays are presented in Sec. 5. Sec. 6 is the summary of our investigations. 



2 Calculations of the conventional mechanisms 

The contribution of the conventional mechanisms are governed by the effective Hamiltonian for 



B decays which is given by [17 



10 \ / 2 10 

neff = k.k; E CrOt + J:c^o, + c,o, + E c^oi + y: c.o. 

L \i=l 2=3 / \i=l i=3 



(4) 



with the effective operators given by 



01 = 


{qaUa)v-A " {upbp)v-A, 


01 


= iMaUp))v-A ■ {Upba))v-A, 


01 = 


{qaCa)v-A " {Cfihfi))v-A, 


01 


= {'laCf3)v-A ■ {Cf3ba)v-A, 


03 = 


{,qJ)a)v-A ■ Y.q'{q'i3(l'(3)v-A, 


O4 


= {qabf3))v~A-J2q'{q'i3(l'a)v-A, 


05 = 


{qaba)v-A ■ Y.q'{Q.'l3Q'p)v+A, 


Oe 


= {qabf3)v-A-J2q'{Q'l3(l'a)v+A, 


Or = 


\{(laba)v-A ■ J2q' ^q' {.^'p'i'p)v +A, 


Og 


= l{Qab/3)v-A ■ T,q' (iq'{Q'i3l'a)v+A, 


O9 = 


\{,qaba)v-A ■ J2q' ^q'{(l'l3(l'/3)v-A, 


Oio 


= I{qabf3)v-A ■ Eg' ^q'Wpl'a )v-A-, 


Og = 


{gj8n') m, q^ a^^ R G% . 






Here q = 


d,s and g'e{u, s, c, 6}, a and (3 are 


the SU{?>) color indices and A^^, A = 1 



(5) 



are the Gell-Mann matrices, and G"^^ denotes the gluonic field strength tensor. The Wilson 



coefficients evaluated at fx = rrib scale are [17 



Ci = 


1.082, 


C2 = 


-0.185, 


C3 = 


0.014, 


C,= 


-0.035, 




0.009, 


Ce = 


-0.041, 




-0.002/137, 


Cs = 


0.054/137, 




-1.292/137, 


Cio = 


= 0.262/137, 




-0.143. 







(6) 



The non-factorizable contributions to i? — > M1M2 can be calculated through the diagrams 
in Fig.l. The results of our calculations are summarized compactly by the following equations 



= % E V;gVja'l{BM^, M2){qu)v-A ® {ub)v-A + a!(5Mi, M2){uu)v-a ® {qb)v- 

V ^ p=U r- ^ 
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+ (5Mi, M2){q'q')v-A ® {qb)v-A + al{BMu M2){qq)v-A ® {q'b)v-A 

+ al{BM,, M2){q'q')v+A ® {qh)v-A + al{BM,, M2){-2){qq')s+P ® {q'b)s-P 

+ a^j{BM^,M2)^eg,{q'q')v+A ® {qb)v-A 



(7) 



+ (-2) ( aliBAh, M2)-e,, + asa{BMi, M^) ] {qq')s+p ® {q'b)s-P 



+ al{BMuM2)h,>{q'q')v_A ® {qb)v-A 



+ [a%{BM,, M2)-e,, + a%,{BMi, M2) j {qq')v-A ® {q'b)v-A 

where the symbol ® denotes (M1M2IJ2 ® ji\B) = (M2|j2|0)(Mi|ji|5). Mi represents the 
meson which picks up the spectator quark through this paper. For Mi a hght vector meson 
and M2 a hght pseudoscalar meson, the effective af s which contain next-to-leading order(NLO) 
coefficients and 0{as) hard scattering corrections are found to be 



al2{BV,P) 
aUBV,P) 

al^{BV,P) 

a3{BV,P) 

<{BV,P) 



a5{BV,P) 
al{BV,P) 



a-r{BV,P) 
as{BV,P) 

asa{BV,P) 
a^{BV,P) 



0, a1{BV, P) = a^{BV, P) = ai{BV, P), i = 3, 5, 7, 8, 9, 10, 8a, 10a. 



C2 


as Cf 


N 


An N 


Ci 


Ois Cf 


N 


An N 


C4 


as Cf 
An N 


N 


C3 


ots Cf 


N 


4^1V 



C^{Fp + Gp{s^) + Gp{s,)) + CiGp{sp) 



+{CA + C^)Y.Cp{sf)+CgCp^g 

f=u 
C5 Cts Cp 



C^G'pis,) + CiC'pis^) + G'pis,)) 



+{C, + C,)j:G'p{sf) + C,G'F^^ 

f=u 



Cg as Cf „ 



12), 



as Cf 
A^W 



{Cs + Cio) E lefC'pisj) + C,^{e,G'p{s,) 

f=u 



ebG'p{sb)) 



i\ Air i\ 



Ots Cp 



3 ^ 3 

f=u 



(8) 



where p = (n, c), g = (d, s), g' = (u, d, s), and / = (u, d, s, c, h). Cp = {N"^ — 1)/{2N) and 
= 3 is the number of colors. The internal quark mass in the penguin diagrams enters as 
Sf = m'j/'ml. X = 1 — X and u = 1 — u. 



Fp 

f. 

Gp,g 

Gp{s,) 



-12\n-^ -18 + fp + fj/, 
nib 



dxg{x)(j)p{x) 



g{x 



3 ma; — 3z7r, 



II 
p 



4vr2 fyf^ 



B 



dz- 



z Jo 



1 — X 

1 



, <Pp{y) 

dy , 

X Jo y 



uux — ie\ 



G' 



P,9 



G'p{s,) 



N A^{0)Ml Jo 

rl 2 

— / dx—4>p(x), 
Jo X 

In 1-4 / dx(j)p(x) / du nnln[s 

3 3 mft JO Jo 

~ lo ^^^'^P^^^ = 

In 1" 3 / dx(j)%(x) / du uu\n [s„ — uux — ie 

3 mi. Jo Jo 



3 rrih 

For Ml is pseudoscalar and M2 is vector, the co-effi cents are 



(9) 

(10) 

(11) 
(12) 

(13) 
(14) 



aUBP,V) 
<iBP,V) 

a^2iBP,V) 

as{BP,V) 

<{BP, V) 



0, a^iiBP, V) = a'liBP, V) = ai{BP, V),i = 3,5, 6, 7, 8, 9, 10, 10a. 



Gs + 



G2 


as Gp 


N 


Att N 


Gi 


as Gp 


N 


An N 


C4 


as Gp 
An N 


N 


G-s 


as Gp 


N 


An N 



GiFy, 

G^Fv, 



Gs{Fv + Gvisg) + Gvisb)) + G^Gyis^) 



+ {GA + Ge)Y.Gv{sf) + GgGv,g 

f=u 



a,iBP,V) 
ae{BP,V) 
a7{BP,V) 



12), 



Ge as Gp 



as{BP,V) 
a,{BP,V) 
aio{BP,V) 

aioaiBP,V) 



as Cp 



07 




iV' 




Cio 




' 


An N 




a, Cf 

l_ 




47r iV 



3 3 



ebGv{sb)) 



where 



Fy = _i21n^-18 + /^ + /^^ 

nib 

fv = dxg{x)(pv{x), g{x) 
Jo 



1 - 2x, 
3 m X — 3z7r, 



II 

V 



'0 

47r2 



1 — X 



N /f-^(0)M| 7o 

/•I 2 
- / dx-(j)v{x), 
Jo X 



dz- 



z Jo 



dx- 



^p[x) 



X Jo 



'dy^. 

y 



2 4 U /"l 

In 1-4 / dx(j)v{x) / du uuln[s„ — 

3 3 nib Jo Jo 



uux — ie\ 



For both M\ and M2 are pseudoscalars , the co-efficents are 

al2{BMi,M2) = 0, a;=(SMi,M2) = <(5Mi,M2),i = 3,5,7,8,9,10,8a,10a. 
<(EMi,M2) 

a^(SMi,M2) 

a3(5Mi,M2) 



<(5Mi,M2) 



a^{BMi,M2) 
al(BMuM2) 



C2 as Cf 
iV+4^1V^^^^- 

^ C3 Cf 

h 

+ (C4 + Ce) ^ C]\i^{Sf) + CgCM2,g 
f=u 

Cq asCp 

T r 

^ iV 47r AT 
^ I , «s Cf 
iV ^i^^AT 



C^{Fm2 + GM2{sg) + CA/2(sb)) + CiGjv/2(% 



■Cq{~Fm2 — 12), 
CiC'p(.,) + C3(G' (s,) + G' (.,)) 



a7{BMi,M2) 



+(C4 + C6)EGm.(^/) + C,C;,2,, 
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as{BMu M2) 
asa{BMi, M2) 

ag{BMuM2) 
aio(5Mi,M2) 

aioa{BMi,M2) 

where 

fm 

f. 



Cs + 



9i 

N'- 



C9 + 

UsCf 



f=u 

GiqFm2) 
G9FM2, 

3 ^ 3 

(Cg + ^'10)2 E ^fGhhiSf) + G<;^-{eqGM2{Sq) + CfeGAfa (■^ft) ) 
f=u 



Clo 


+ 


as Gf 






G9 


+ 


as Gf 


N 


Att N 



_121nA_i8 + /£ 
1 



M2 





47r2 



ruh 

dxg{x)(t>M2{x) 
fh'h fs 



Ah 



rll 
J M2I 



GM2,g 
Gm2 {Sq) 



N f'^^^^'{Q)Ml Jo 

- / dx-(j)M2{x) 
Jo X 



dz- 



9{x, 
z 



= 3 111 X — 327r, 

1 — X 

'dx^^ tdy'^. 
X Jo y 



G'. 



3 3 m;, 
"1 



M2,g 



G'u^iSq) 



1 

3 



dx7^4>M2{x 
In A + 3 



4 / dx(t)M2 (3^) / "^"^ ''^'^ 111 [■So — uux — ie] 
Jo 

_3 



uux — ie\ 



(21) 

(22) 

(23) 
(24) 
(25) 
(26) 

mfe Jo ^0 

where (f)p{x) and (f)^p{x) are the pseudoscalar meson's twist-2 and twist-3 distribution am- 
phtudes (DA) respectively. 0y(x) = (j)v,\\{x) is the leading twist DA for the longitudinally 
polarized vector meson states. We have used the fact that light vector meson is longitudinally 
polarized in i? — PV decays. In the derivation of the effective coefficients Oj's we have used 
NDR scheme and assumption of asymptotic DAs. The infrared divergences in Fig.l.a — d are 
cancelled in their sum. 

With the effective coefficients in Eqs.^, |TB| and we can write down the decay amplitudes 
of the decay modes (we only list four decay modes here which have been observed by CLEO. 
The other decay modes are given in appendix A) 



1 nl 

dx(j)\.j^ (x) / du uu In [s^ 

•J 



mlWcbil - Rce'^ ai{X) + <(X) 
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+RkK{X) + a^{X) + a^a{X)] + aio(X) + aioa{X) 
+al{X) + aio(X) + aioa(X) + RK[al{X) + aglX) + Ogal^)] 



-^a8(r) + a8.(y))(/,^, - + \{-aj{Y) + a9(F))(/,", - J^,) 

+ (-iaio(F) + aioa(V'))/,^l + (o3(r) - a5(r))(2/;, + f^,) + al{Y)f^ 



+i?^,(a^(r) - -as{Y) + a8,(F))(/,^, - /;,^ 



+ -(a,{Y) + ar{Y)){f;,-i:;,] 



-a,o{Y)-awa{Y))K 



(27) 



mt/i X = B-r]',K- and Y = B-R-,?]'. 



+ i(-a7(X) + a9(X))(/;, - /^^) + (-^aio(X) + aioa(X))/,^ 

a3(X)(2/;, + f^,) + al{X)f^, - a5(X)(2/;, + f^,) 
+ Rt,.{al{X) - ]^a,{X) + aUX)){r,, - f^,) 
+ i(a9(X) - a7(X)+)(/^^ - /;,) + (-iaio(X) + a,oa{X))f^, 

I • / /■ 2 

+ Rxia'^iY) - ^as{Y) + asa{Y)) - ^aio(F) + a,oa{Y) 

+ al{Y) + RKKiy) - l^^siY) + asa{Y)) - ^a^oiY) + a.oaiY)^ 
with X = B^K^, r] and Y = 5%', 



(28) 



M{B^^7^K*+) = %<-^*(mJ)m|V;,(l-^)^i?,e^^ a2(X)/; + a3(X)(2/; + /- 



v^^ 



+ a:{X)f; - a5(X)(2/; + f;) - i?;«(X) - -a8(X) + a8.(X))(/^^ - /;) 



+ -(a9(X)-ar(X))(/;-/- 



^aio(X) + a,oa{X))f^ + a3(X)(2/; + f^) 
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+ i(a9(X) - a7(X))(/; - f^) - (iaio(X) + aioa W)/,^} 

+ aio(F) + aioa(l^)]+a4(>^) + aio(n + aioa(>^)}, (29) 
X = B-K*-,r]' and Y = B-r]',K*-. 



M{B° r]K 



*0\ 



+ 



+ al{X)f:; - a,{X){2^ + r,) - W^{al{X) - -a,{X) + aUXW, - f^) 



-(a9(X)-a7(X))(/;-/- 



:-aio(X) + aioa(X))/, 



a3(X)(2/; + /- 



+ -{a,{X) - aj{X)){f; - p - (-aio(X) + aioa(X))/,^ 



<(y)--aio(l^) + OiOa(>^) 



a^(>^)-^aio(n + aiOa(n}, 

mt/i X = B°K\7]' and Y = B^r]',K*^. 



(30) 



Where Rr 



of v:,. X 

in Sec. 4. 



IK,. 



A I y^^^ g^j^^ y^^^ g^j,g chosen to be real and 7 is the phase 



1-A2/2 I Kb 

0.2196. We will present inputs and numerical results for those magnitudes 



3 Re-examination of the spectator-hard-scattering mech- 



anism 



3.1 Calculation of g*g* — t]^'^ transition form factor Fg*g*_^{o{Qi, Ql) 

Recently Ref. |T9[ studied g*g — r]^'^ transition form factor Fg,g_^(,){Ql,Ql = 0). The same 
method can be easily used to calculate Fg,g,_^{i) {Q\, Q^), where both of the two gluons are off- 
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shell. To keep the completeness of this paper we recapitulate the main points of this calculation 
here. 

The 77' meson is mainly a flavor singlet meson. Because of its singlet structure the 77' meson 
may have gluonic content through the QCD anomaly. The contribution of gluonic wave function 



to the transition form factor Fg,g_^(i){Q\,Q\ = 0) has been tested, which is very small |jT9|. 
By analysing the gluonic wave function, when both of the two gluons are off-shell, it will give 
an extra scale suppression. So for the case of Fg*g,_^{i){Q\,Ql), it may be safe to neglect the 
contribution of gluonic wave function of 77' meson. 

We take the qq — ss mixing scheme for the t]*-'-* meson in this calculation. Here qq means 
the light quark pair uu and dc? |Q . In this mixing scheme the parton Fock state decomposition 
can be expressed as 

I rf) = cos(f> I rig) — sincj) \ ris), 

I T]') = sincj) I ?7g) + coscj) \ rjs), (31) 

where is the mixing angle, and | ?7g) ~ | + dd), \ rjs) ~| ss). 

The diagrams contributing to the transition form factor (Qi, Qi) shown in 

Fig. 3. By direct calculation of these two diagrams the g*g* — 77*^'-' transition form factors can 
be parameterized as 

< glQllv^''^ >= -4na,6abte^,^pQ'^Q^2Fg.g,_,io{Ql, Ql). (32) 
and Fg*g*_^(,){Ql,Q2) is found to be 

Fg*g*_r,('){Ql,Q2) = ^ f^C) '^^7^ 



q=u,d,! 



x)Qi + xQl — x{l — x)m^(,) + ie 

(33) 

where the variables /^(,) can be related to the decay constants of | rjg) and | rjg) 

f^> = ^ sin 0, f^, = fs cos 0, 

/; = ^cos0, /,^ = -/,sin0. (34) 

The decay constants fg, fs and the mixing angle have been constrained from the available 
experimental data, fg = (1.07 ± 0.02)/^, /, = (1.34 ± 0.06)/^, = 39.3° ± 1.0° [§. 
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To the accuracy of this paper, (x) is taken to be the leading twist DAs as [|T^ 0^(,) (x) = 
6x(l — x). 

The transition form factor will play a pivotal role in estimations of gluonic exclusive pro- 
duction of T]' and rj. 

3.2 The magnitude of spectator-hard scattering mechanism for B 

rj'M 

In this subsection, we would re-calculate the magnitude of spectator hard scattering mecha- 
nism(SHSM ) for i? ^ r]'K proposed in Ref. 0] and generalize it to -B ^ rj^'^M where M is a 
light pseudoscalar or vector meson. 

The SHSM is described by the Feynman diagrams in Fig. 3 where b quark decays to s(d) 
quark and a hard gluon. Since the virtuality of the gluon is much larger than the typical scale 
of QCD boundstate Aqcd i-^., l/Qf << I/Aq^^, it would fluctuate into a small size fast flying 
color-octet qq pair well before it fiys out of the light cloud of B, so it would hard interact with 
the spectator. More arguments for validity of using perturbative QCD in this calculation could 
be found in Ref.0. 

For M a light pseudoscalar, the amplitudes Mi for Fig.3.a and M2 for Fig.3.b are calculated 
to be 

Ml = -^{-V:^V,,)Cgm,fBfp{2f;,o + r,vMl^ x (m^Fi^a - /iP^t^^s), (35) 



Ftl = -8/ dzdyMz)rp{y) 



dx—n, — o -• (36) 

xQi + XQ2 — xxm^{,) + te 



^2 



f£1 = -8 / dzdyM^)<Pp{y) 
Jo 



— -V-^r>'V-/^_,. ^g2(g2 



M2 



dX _ r, „r, _ o . , ("^''') 

xQi + XQ2 — xxm^f,) -|- le 

^ E v;,Vf,c(fBfp{2f;,, + /,^))«.^^ x f^;; (.^), (38) 



Kpisf) = C dzdyMz)rp{y)^[Q\-^^''^ 
Jo y - z \ 



2) 



y- z \ Q 
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/ duuu \ —1 + 2\n. ln(sf — MM — \ — ie) 

Jo \ rrih mr 



1 </)^f,) (x) 

xQi + x(y2 — xxm^(i) + ^e 

When M is a light vector meson, the amphtudes M3 and M4 for Fig. 3. a and Fig.S.b respectively 
are 

<Vi = -8 I'dzdyM^Wiy)^-^ i n?(n2^f^ - ^] 
Jo y- z V<5i(<52 + «e) / 



1 (p^'ix) 
xQi + xQl — xxrn?^,) + ie 



X / dx——K -75 — — — K -. (41) 



M4 



^ E V;,Vj,C(Mvi2f;,, + /,^,)a^^ x F^p,.2(^/), (42) 

^^pV2(^/) = 16 f'dz fdyM^Wviy)^ [qI - ^^'-^'^' ] 
Jo Jo y — z \ V2 / 

X / dnnn — l+21n \n(sf — uu — \ — ie) ] 

Jo \ nib m( J 

X / dx—n^ — ^ -. (43) 

JO xQi + XQ2 — xxm^(,-j + le 

Where q = d,s and f = u,c. Qi = {1 — z)Pb — (1 — y)PM and Q2 = zPb — yPu- 

It should be noted that the effective vertex h sg* calculated from full theory in Ref. 
could not be used here, otherwise the contribution of top quark penguin will be double counted 
when amplitudes of SHSM are added to these in Eq.^. Other four fermion operators, say, 
O3, O4 and Oq can also be inserted in Fig.3. However, dominate contributions come from the 
insertion of O"''^. 
With note of 



and the DAs to be given in next section, we can see that Fj^^l, Fpp (sj), -Fepyi, and F'ppy2{sf) 
are free of IR divergence. Since the twist-3 DA (j)%{y) for pseudoscalar meson does not fall 
at end point, integration oil/{y — z) over y results in ln(2), but such large logarithm will be 
smeared by DA of B meson through the integration / (psiz) In zdz, which are shown by analytic 
expressions in Appendix. B. Generally the integrations are free of IR divergence. 
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4 Numerical calculations and discussions of results 



In the numerical calculations we use \Wl 



t{B+) = 1.65 X lO^'^s, r{B") = 1.56 x 10~'^s, 

Mb = 5.2792^61^, mf, = 4.8GeV, = lAGeV, 

rris = SOMeV, fs = O.lSOGeV, U = 0.133Ge1^, 
Jk = 0.158GeV, fx* = 0.2UGeV, fp = 0.21GeV. 
For the chiral enhancement factors for the pseudoscalar mesons, we take 

i?7r± = Rk±'0 — 1-2, 



which are consistent with the values used in |14, 22, 23 1, and 



m 



,(') 



{') 



msTnii 2 

We take the leading-twist DA and the twist-3 DA (f)^{x) of light pseudoscalar and 
vector mesons as the asymptotic form |^ 

(44) 



0py(a;) = 6x(l — x), </'p(x) = 1. 



For the B meson, its DA is modeled as [25, 2f: 



^b{x) = Nbx'^{1 — a;)^exp 



2^1 



(45) 



with ujb = 0.3 GeV, and Nb is the normalization constant to make Jq dx(f)B{x) = 1. (f)B{x) is 
strongly peaked around x = 0.075, which is consistent with the observation of Heavy Quark 
Effective Theory that the wave function should be peaked around Aqcd/Mb- We get the object 
Jq dx(f)B{x)/x = 14.62 which is very near to the argument of Jq dx(j)B{x)/x mB/ \b = 17.6 in 
Ref.0. 

We have used the unitarity of the CKM matrix V*gVub + V*qVcb + V^*Vtb = to decompose 
the amplitudes into terms containing V*^Vub and V*^Vcb, and 



|\4^| = A = 0.2196, iKb/Kbl = 0.085 ± 0.02, 
|K.| = 1-AV2 |V;d| = l-AV2. 



(46) 
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We would use the latest CLEO results for |Vcb|[g3, EH 

\Vcb\ = 0.0464 ± 0.0020(stot.) ± 0m21{syst.) ± 0.0021(t/ieor.), (47) 

and leave the CKM angle 7 as a free parameter. 

For the form factors for B n, K, K* and p, we take the results of light-cone sum rule 



29, 30 



F^-"*(0) = 0.3, F^^^(O) = 0.36, A^"'^ (0) = 0.372, A^^^*(0) = 0.470. (48) 

Compared with these rather well studied form factors, the form factors for B rj^'^ are poorly 
known, which has hindered theoretical predictions for B decays involving r/^'^ very much for 
a long time. Neglecting 77*^'^ masses effects, we argue the following scaling relations for from 
factors at large recoil point in the heavy quark limit nib 00 

F^-'-{0) A' F^--(O) // ^ ^ 

We get 

^B-r,^0) = 0.176, F^^"'{0) = 0.142. (50) 

which agree well with the values F^^''(0) = 0.181 and F^^'''(0) = 0.148 in Ref.[||, ||. 

Taking 7 = 55° as benchmark, we present our numerical results in Table. 1. As shown in the 
table, for the case of F^'^^'^O) = 0.176(0.142), our estimations of the conventional mechanism 
contributions to B^ K'^t]' and B^ — * rj'K^ are about (30 ~ 40) x 10"^ which confirm many 



theoretical estimations in the literature p2|, It implies that significant new contributions 



to those decays are needed to interpret the CLEO data. Before going further, we note that 
the theoretical predictions for B rjK* are also much smaller than their experimental results. 
In contrast to i? — > r]'K, new contributions to B decays rjK* from both b s{cc) si] and 
SHSM are negligible. Because 77 and rj' mesons are much heavier than vr, the recoiling of r]^'^ 
should be smaller than that of vr, so the scaling relation of Eq.49 might be broken to certain 
extent. Driven by the facts, we pose 30% enhancement to the form factor F^"*'^'^'\o). From 
column 4 of Table. I, we can see the 30% enhancement of p^^'n^'^ (Q) is preferred by experimental 
data. However, conventional mechanism contributions to B^ r]'K^ and B'^ rj'K^ are just 
enhanced by 17% and still around 1/2 of their experimental center values. With incorporation 
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Table 1: Branching ratios (in units of 10~^) for B charmless decays involving 77^'). Experimental 
results are taken from |TT|. Our results are made for four cases with 7 = 55° and two cases with 
7 = 85°. The labels CM A and CM B represent estimations of the conventional mechanism 
contributions(CM) with inputs F^'^^'\0) = 0.142(0. 176)(A) and F^''*''(0) = 0.185(0.229)(B) 
respectively. Columns with CM A(B)+SH represent our results with the incorporation of 
spectator hard scattering mechanism. 



Decay 




7 = 55° 


7 = 85° 


CLE0[11 


1 


Modes 


CM A 


CM A+SH 


CM B 


CM B+SH 


CM B 


CM B+SH 


B or 90% UL 


B+ 


ri'K+ 


27.1 


63.4 


32.8 


68.9 


35.4 


69.3 


801.^° ± 7 


T->0 




28.4 


67.1 


35.0 


74.8 


34.5 


73.1 


89tll ± 9 


B+ ^ 


r]K*+ 


5.40 


5.24 


6.88 


6.73 


9.62 


9.44 


26.4l8:2 ± 3.3 




r]K*^ 


7.40 


7.27 


9.96 


9.79 


9.82 


9.65 


IS.Stli ± 1-6 


B+ ^ 




3.78 


7.60 


6.12 


10.8 


4.86 


11.0 


<12 


50 ^ 




0.13 


0.51 


0.16 


0.52 


0.20 


0.74 


<5.7 


B+ 


7]'K*+ 


5.87 


6.40 


5.02 


5.47 


3.92 


4.21 


<35 


B^^ 




3.22 


3.68 


2.10 


2.49 


2.06 


2.42 


<24 


B+ ^ 


ri'p+ 


3.73 


3.71 


6.36 


6.33 


6.54 


6.60 


<33 






0.020 


0.010 


0.023 


0.012 


0.033 


0.020 


<12 


B+ ^ 


rjK^ 


3.11 


7.85 


1.49 


6.00 


1.21 


4.69 


<6.9 




r]K^ 


1.64 


8.89 


0.28 


8.88 


0.28 


7.87 


<9.3 


B+ ^ 


?77r+ 


6.54 


4.84 


10.2 


8.02 


7.96 


5.15 


<5.7 


50 ^ 


rjTT^ 


0.36 


0.41 


0.42 


0.46 


0.53 


0.60 


<2.9 


B+ ^ 


rip+ 


6.45 


6.39 


10.7 


10.6 


10.3 


10.4 


<15 






0.02 


0.017 


0.023 


0.001 


0.033 


0.015 


<10 
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of SHSM contributions, theoretical results are improved much, and the experimental data from 
CLEO could be accommodated, although our results are still slightly lower than the center 



values. It is also noted that BaBar collaborations |]32[ have reported their preliminary results 

B{B+ T]'K+) = (62 ± 18 ± 8) X 10^^ B{B^ r]'K^) < 112 x 10"^ (51) 

where the center value of B{B^ rj'K^) is smaller than that of CLEO's. It might be safe to 
conclude that the large branching ratio oi B ^ rj'K could be understood in the SM. 

Topologically, SHSM would contribute to the decays B rj'K* the same as it to -B ^ r]'K. 
However, SHSM amplitudes depend on the spin configurations of K and K*. As shown in 
i?gs.|3^ — compared with the SHSM amplitudes for B K*ri', SHSM amplitudes for 
B rj'K are chirally enhanced by K twist-3 DA. So that the new contributions to i? ^ rj'K 
are much larger than its to 5 — > rj'K*. In Fig. 4, we display our results for the branching ratios 
of B two-body charmless decays involving r]^''> as functions of the weak angle 7. From Fig. 4, 
we can see that SHSM contributions to B rj' + P are much larger than its contributions 
to B 7]' + V. It is easy to understand that SHSM contributions to B rj + M are very 
small because of cancellations between (uu + dd) and ss contents of rj. For the four measured 
decays B~^ {rj'K~^, rjK*) and B^ {ri'K^,riK^*), our results generally agree with CLEO 
data. In our calculations, we have used large form factor F^^^{0) to give large B{B^ rjK*^) 
and B{B^ rjK*^), meanwhile we meet constraint from 5+ 777? whose upper limit set by 
CLEO is B{B^ rjTx^) < 5.7 x 10~^. Since B^ rjn'^ is dominated by the tree amplitude 
with z/^F'^^''(m^)ai, B{B^ rj-K^) is very sensitive to the form factor F^^^(m^). To meet 
the large branching ratio B{B^ rjK*'^) = {26.4:1il'^ ± 3.3) x 10^^ and comply with the small 
upper limit B{B^ rj-K^) < 5.7 x 10~®, we need cos 7 < to enhance B{B^ t]K*^) and 
suppress B{B^ rj-K'^). This situation is very similar to that in B^{B^) vr^vr''^, K^n^ and 
B 7rp,7TK* decays as discussed in Refs. [p^ , |33| , |3^ , |36| , |37[| recently. It is worth to note 



that most recent theoretical analyses of CLEO data based on different approaches endorse, 
although not definitely, negative cos 7 to some extent. However global CKM fit has given the 



constraint 7 < 90° at 99.6% C.L. [^, With refined measurements at running B factories 
BELLE and BaBar to come very soon, the following B exclusive decay modes will give strong 
constraints on 7 

PP modes : 5° ^ tt^tt^, 5° ^ A'^vr^, B^ K^n^, B^ r/V^, ^ ^7r±; (52) 
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PV modes : 5° vr^p^, n^oo, 7r±p°, B^ K^oo, B^ r]K 



*± 



(53) 



Branching ratios for the above decay modes are of order of 10"^ ~ 10^^ which can be well 
studied at B factories to constrain cos 7 tightly. If the disagreement of constraints on 7 between 
global fit and direct model calculations becomes serious, it might imply the failure of the models 



employed here and in Refs. p2|, |33| , ^ |3^. Very probably theories for B hadronic decays 
will be tested and driven by the observations to be made at BaBar and Belle. 



5 Final states interactions and CP violation 



As shown by Beneke, Buchalla, Neubert and Sachrajda in Refs. |p!4|, |I5|, the final states in- 
teractions in charmless B two-body decays are calculable in the QCD improved factorization 
framework, which turn out to be nonfactorizable corrections. The nonfactorizable corrections 
for the decays studied in this paper are shown in Eqs.^, SHSM amplitudes are generally 

nonfactorizable and always contribute large strong phase to certain decay modes. It is worth 
to note that some shortcomings in the "generalized factorization" are resolved in the frame- 
work employed in this paper. Nonfactorizable effects are calculated in a rigorous way here 
instead of being parameterized by effective color number. Since the hard scattering kernels are 
convoluted with the light cone DAs of the mesons, gluon virtuality k'^ = xml in the penguin 
diagram Fig.l.e has well defined meaning and leaves no ambiguity as to the value of k'^, which 
has conventionally been treated as a free phenomenological parameter in the estimations of the 
strong phase generated through the Bander, Silverman and Soni(BSS) mechanism ||40||. So that 
CP asymmetries are predicted soundly in this paper. 
The direct CP asymmetry parameter is defined as 

CP \M{B+ ^ f)\^ + \M{B- ^ f)r ^ ^ 

and 

CP |_M(EO^/)|2 + |A^(fio^/)|2- ^ ^ 

For CP-violations in B'^{B^) 'q^'^^ir'^^'q^'^p^, B^ — B^ mixing effects should be considered. 
However, the branching ratio for these decay modes are very small( below 10^^). We would 
not address CP violations in those decay modes in this paper. 
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Our numerical results for the direct CP violations are shown in Fig. 5 as functions of 7. For 
7 e [50°, 90°], the direct CP-violation in rj'K^ are found about 8% ~ 10%. The recent 

search for direct CP violation in — > rj'K^ made by CLEO Collaborations has reported 



A'^JpiB'^ r]'K^) ~ +0.03 ± 0.12 (56) 

which results in the 90% C.L interval [-0.17, 0.23] for A'^'^'^B^ -> rj'K^). With more and better 
data to come soon, the sensitivity of A'^"'{B^ rj'K^) at B factories could reach ~ ±4%. 
The prospect of observing direct CP violation in B^ r]'K^ are expected to be quite good. 
The direct CP violations in B^ rj'K'^, r]K* and B~^ — s> rjK*'^ are also estimated to be about 
few percent but smaller than A'^^'^{B^ rj'K^). Considering their experimental sensitivities 
and/or branching ratios, prospects of observing direct CP violations in these decay modes are 
much weaker than in B^ — * r]'K^. 

Dighe, Gronau and Rosner|H2[ have predicted large CP violations in B^ rjir^ and B^ — > 



r]'TT . More earlier similar conclusion could be found in Ref.[^. For 7 G [50°, 90°], We find 



A^}>{B+ ^ r/'7r+) ^ -50% ~ -80%, A^Sp{B^ ^ V^^) ~ +15% ~ +30%. (57) 

From Fig.5.5 and Fig.5.11, we can see that SHSM contributions could enhance ^^*^(i?+ — >• 
ry(')7r+) very much. As shown in Fig.4.5, the branching ratio B^^^^B'^ rf'-K^) is predicted to 
be of order of 10~^. 

It should be very promising to observe direct CP violation in B^ — > r/V^, r/vr^ in the 
near future. We also predict large A'^^'^{B'^ t]K^). The decay modes get a large strong 
phase through SHSM as shown in Fig. 4. 11. The strong interference between tree and penguin 
amplitudes leads to 

J^d^r^Q+ _^ ^^+^ _ _^20% ~ +50% (58) 

for 7 e [50°, 90°]. 

To summarize this section, we find large direct CP violations in decays B"^ v' ~^ 
rj'ir^, rjiT^ and B^ r]K^. 
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6 Summary 



With the newly developed QCD improved factorization approach [p!4|, p!5|] , we have studied 
hadronic charmless two-body decays of -B„ and Bd involving rj oi rj'. Nonfactorization effects 
are calculated in terms of order of 0{as) corrections to the leading factorizatable amplitudes. 
We find again that the conventional mechanism account for about one half of the decay rates 
of -B — > rj'K. Significant contributions are needed to solve the "puzzle" of unexpected large 
branching ratios. 

To clarify possible sources of new significant contributions to B{B "rj'K) up to date, we 
display the following facts. Theoretically new contributions due to intrinsic charm content of 
the rj' have been realized to be small. Motivated by their observation of large B decay rates to 
?7'i^[|l], [n]], CLEO Collaborations have made searching for B decays to ?7c/'^||2^. It is found that 



there is no unexpected enhancement in the rjcK rate. Since the rjcK rate should also be enhanced 
if the rj'K rate enhanced by the intrinsic charm content of the 77', the CLEO results indicate 
that the charm content of the rj' may be not the explanation for the anomalous E{B fj'K). 
Alternatively one may turn to new physics. We have known that the experimental observations 
of -B — > and B K*'~f agree with the SM expectations, which implies that there is no 
anomalous large new physics effects in B decays. Motivated by these considerations, we have 
re-examined the SHSM. Starting from calculations of the transition form-factor g*g* — 77', we 
have incorporated the new contributions to that of conventional mechanism for B decays and 
successfully accommodated the large rate of B decays to rj'K in the SM. We have found that 
the absolute strength of the amplitude of SHSM for B decays to rj'K are about 1/3 ~ 1/2 of 
the absolute strength of the conventional mechanism which shows the conventional mechanism 
to be the dominant. As shown in detail, the SHSM is important for B decays to rj'K but not 
for B decays to rj'K*. 

We have estimated the direct CP violations in the decays. rj'K^) is found to 

be around -|-8% ~ -1-10% for 7 G [50°, 90°]. Due to significant contributions from SHSM, we 
predict 

^dir^^± _^ ^/^±) _ 4Q^^ _ 7Q0/^_ A^'^-iB^ ^rin^)- 20% ~ 40%. 

We also predict large direct CP violation in B decays to rjK^. Prospects for observing direct 
CP violations in B^ decays to rj'K^, rj'-iT^, r]7r^ and rjK^ are quite promising at the ruining B 
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factories BaBar and BELLE. 

Note added After we finished this work, we note the work[|4|| by Ah and Parkhomenko 
which is focused on g*g* — t]' transistion form factor. They find that the contribution gluonic 
content of rj' to the form factor Fg*g*rj'{Ql, Ql) as large as few ten percentages for smah and 
rather small for large Q^. Our formfactor agree with theirs to leading terms. Additionally, in the 



framworkp!4|, |T5[ employed here all diagrames in Fig.l and Fig. 3 are subleadling nonfactorizable 
contributions, so the gluonic contributions could be neglected. 
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Appendix A 



The decay amphtudes the conventional mechanism of some of the decays in terms of the 
effective coefficients a^'s are presented as the followings, 



MiB+ -> rj'n+) 



+a^^{B+n+,v')Q,-, + [a^{B+n+,v')-'^as{B+n+,r,') + asa{B+n+,v')]R''^^„Q„ 

+ i[a9(B+^+,V)-a7{B+^+,V)] {/;'(,) " /^w ) + [-^«lo(B+^+, V) + aiOa(B+^+, r,')]/;,,) 

- [as{B+7r+,rj')-a,{B+7r+,v')]{2Q,^+f^,„) + al{B+7r+,r,')Q,^ + [al{B+7r+ ,r,') 

-ia8(S+^+,r,') + a8a(B+^+,»?')]fi^(,)/> + h{B+^+, V) " ^7(6+^+, r;')] ^ (/;<„ " ) 

\i{B+r,',n+) + a^{B+r,',n+) + [a^{B+v', n+) 



+i^(M| -m2(,))FS^'''''{0)/^AV;6<; Rue 
+asiB+ri',7v+) + asa(B+r,', R^ + aw{B+n', 7r+) + awa{B+ ri' , 7r+) 



al(B+ri',^+) 



+ [al{B+r)' + a8(B+V, 7r+) + asa(B+ ri' , 7r+)] R^+aio(B+ V, ^+) + awa{B+ -q' , tv+ ) 



(59) 



M{B° ^ rjCvrO) 



M{B+ rj'K*+) 
M{B° ^ri'K*'^) 

M{B+ ^p+»7<')) 



-^^MiF«— (m2,„)Ay,,|R„e'- [a2(B«.°, r^W )/;<„ + [aaCfi".", ryW) - as (5°.°, r,(') )] (2/;,,, 
+f;,,,)+al(B'\\r,'~'))f;^,, + [al{B^\\rf'))-\a,{B^\\^^^^^ 

+ i [a9(B0.^r,(')) - a7(B«.^r,W)] (/;,,,-/;,)+ [-iaio(B°-^r,(')) + aio.(B«-^r,W^ 
"[a3(B°^°,r,W) - a5(B°^°,ry('))] (2/;(„ +/;<„) + r^O)/;,,, 
+ [ag(iJ«7r°,r,('))- la8(B0^°,r,(')) + a8,(iJ«^«,^W)] fl^,,, 

+ i [a9(B'V^r,(')) - a7(B'V^r,('))] /^,,)) + [-iaio(B'V^^W) + aioa(B°-^r;W)] 



2 ^ 



II 



^ 7?it*+) iv rf); 
M{B+ ^ r]K''^) (r, ^ r,'); 



(60) 

(61) 
(62) 



2(B+p+,r,W)/;,„ + [a3(i?+p+,.?W)-a5(B+p+,r,W)] (2/;,,, +/;(„) 
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M{B+ riK+) 

where i?„ 



+a^(B+p+,r,W)/^"(,) - [<(B+p+,r,('))- ia8(B+p+,r,W) + a8„(B+p+,r,('))]R^„)/^"(,) 
+ l[as(B+p+,ri'-'^)-a,{B+p+,v<-'^)]{Q,,- r^,,,)+[-'^au^ 

- [a3(B+p+,r,('))-a5(B+p+,r,('))](2/^"„) + ) + ^(') ) _ [ag(B+p+ , r,(') ) 

+ [- iaio (B+ p+ , W ) + aioa (B+ P+ , r; W )] 

'ai (B+r,W , p+) + (B+ryM , p+) + aio (B+r)''' , p+) 



+aiOa(B+r;('',p+) 



(B+ri^'') , p+) + aio (B+t^*'' , P+) + aiOa (B+r)''' , p+) 



(63) 



-^MK^^(m2„)AV;,<jiJ„e'- 



a2(B«p«,V'))/-„ + [a3(B0p«,V'>)-a5(B°P°,r,('))](2/;(„ +/;„) 

+aUB°p'>,v^'^)Qo - [a6(i?%°,'?"')- ia8(B»pO,r,(')) + a8.(B«p«,„W)]B^„/;„) 
+ 1 [a9(BV,r;(") - a7(B0p°,^(")] (/;„) - /;,„)+ [-iaio(B«p°, r^O) + aio.(B%''/^W)] /^j,," 
-[(a3(BOp°,r,W)-a5(B0pO,r,('))] (2/;(„ +/;„) + aS(B«pO, r^W )/;„ - [a^(B«p«, r,W) 
-ia8(B0p«,„(')) + a8jB0p«,„(")] + [a,{B° p^V^'^) - ar(B«pO, „('))] (/;,, ~ /;„) 

+ [-iaio(B«p«,„(')) + aioa(BOp«,r, ('))]/;,„ 



+ ^Mlf,F^-'^"' (0)AK, <j L(B%('), p«) - a,"(B%(') , p«) + | [a9(B%(') , p") + ar(B%W, p^)] 



-iaio(B^('\p«)-aioa(B%(",P») 



-iaio(BV',P°)-aiOa(B°r,W,P°) 



-aS(B%(') , p«) + I [ag (B°r,(') , p") + a7(B%(') , p«)] 



M{B+ -> r;'_ft'+) (V j?); 
A4(B+ ^ r?'i^O) (V ^ »)). 



(64) 

(65) 
(66) 



M2 



(') 



(1 



%(') 



and the corrections to the decays from 



SHSM can be read from equations in Sec. 3. 2. 

Appendix B 



and Fgpyi 


are g 


F^"' - 

^tw2 


Jo 




Jo 


^ BPVl ~ 





dz 



z 

b{z] 



z{l-z) 



Yi{z, /c^(,),m^(,)), 
Y2{z,kj^(,),m^(,)), 



(67) 
(68) 
(69) 
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and k^(o = m\/m?^(,y The functions Yi{z, k,m) and Y2{z, k,m) are given by 

1 



Yi{z, k, m) 



Y2{z, k,m) 



[k — 1)^ m? 



2k (1- z) - k^ {1- zf + 2 (l-kz) 
+3k {1- z) il-kz) -2{l-kzf -2 ln(l - k {I - z)) 
+3 A; (1 - z) ln(l - k {I - z)) - P (1 - zf ln(l - k {I - z)) 
+k^ (1 - zf ln(-A; (1 - z)) + 2 ln(A; z)-?,k {I- z) ln(A; z) 
+3 A; (1 - 2) {l-k zf hi{k z)-2{l-k zf ln(fc z) 
-3k (l-z) {1-kzf ln(-l + A;^) + 2(l-A;^)^ ln(-l + A;2) 
1 

{k — 1) 

- ln(l - k {1 - z)) -2 {1 - z) ln(l - k {1 - z)) 
+2 A; (1 - z) ln(l - k {1 - z)) + 2 k {1 - z f ln(l - A; (1 - 2)) 
-A;2 (1 - zf ln(l - A; (1 - ^)) - 2 A; (1 - ln(-A; (1 - ^)) 
+k^ (1 - ln(-A; (1 - ^)) + 2 (1 - z) ln(l - k {1 - z) - z) 



3 -1 + A; {l-z) + kz 



-2 A; (1 - ;z) ln(l - k {1 - z) - z) 



-2{l- zf ln(-A; (1 - ^)) ln( 



l-k{l-z)-z. 



+2k{l-zf \n{-k {l-z)) ln( 



I- z 
l-k{l-z)-z^ 
1 - z ' 



+2 (1 - zf ln(l - A; (1 - z)) ln(-^^-^i^^^) 
-2 A; (1 - zf ln(l - A; (1 - z)) ln( ~^^~^ (1 " ^) " ^) ^ ^ j^^^ 
+2 (1 - z) ln(A; z) - 2 A; (1 - 2) ln(A; z) - 2 {I - z) {I - k z) ln(A; ;2) 
+2 A; (1 - z) (1 - A; z) ln(A; z) - (1 - A; z)^ ln(A; z) 
+2 (1 - ^) (1 - A; ^) ln(-l + A; ^) - 2 A: (1 - ^) (1 - A; ^) ln(-l + kz) 
+ {l-kzf ln(-l + A;^)-2 {l-z) \n{-z + kz) 
+2k (l-z) \n{-z + kz) + 2{l-zf ln(-l + A;z) ln(-~^^^^^ 



-2A;(l-z)Mn(-l + A;z) ln(- ^ + ^^^ 



+2A;(1 - zf \n{kz) ln(- 



1 - z 
-z + k z) . 



1 - z 



2{l-zf ln{kz) ln(- 



(70) 



-z + k z) 



zfU2{k) 
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+4 (1 - zf Li2{k) + 2 (1 - zf Li2 I 
-2Hl-zfa.iZS-l±Jl^^- 



-2k{l-zfLi, 



(-l + Ml-z)), 



+ 2{l-zfLi2 



'l-kz' 
l-z 



(71) 



Where Li2{x) is the dilogarithm or Spence function. 
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Figure captions 



FIG.l Order as non-factorizable contributions in 5 — M1M2 decays. 

FIG. 2 Tlie triangle Feynman diagrams for the transition form-factors of g*g*—ri' and g*g*—ri. 

FIG. 3 The Feynman diagrams of the spectator hard scattering mechanism for B decays to 
r)'M, where M is a light pseudoscalar or vector meson. 

FIG. 4 Branching ratios for B decays are shown as curves as a function of 7 in units of 
10^^. The dashing curves are the results of conventional mechanisms, the solid curves are our 
estimations with incorporating SHSM contributions. The branching ratios measured by CLEO 
Collaboration are shown by horizontal solid lines. The thicker solid horizontal lines are its 
center values, thin horizonal lines are its error bars. 

FIG. 5 Direct CP violations in B decays involving ri^'\ Dashing curves are the predictions 
of conventional mechanism estimated with the QCD improved factorization approach. Solid 
curves are results when SHSM contributions added. 
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Fig.4.1, Br{B+ ri'K+) vs. 7 
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Fig.4.4, _Br(_BO -> jjE"*") vs. 7 
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Fig.4.13, Br{B+ ^rin+)vs.^ 
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Fig.4.2, Br{B° r/' K°) vs.'y 
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Fig.4.5, Br{B+ r^'-K+)vs.^ 
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Fig.4.8, Brl^B^ -» rf K*^)vs.'y 
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Fig.4.11, Br(B+ -» r]K^) vs. 7 
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Fig.4.14, Br(BO Tjn°) vs. 7 
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Fig.4.3, Br{B+ riK*+) vs. 7 
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Fig.4.6, Br{B^ ri'Tr'^)vs.^ 
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Fig.4.9, Br{B+ ~* rj' p+) vs.'y 



50 100 150 200 250 300 350 

Fig.4.12, Br{B° riK°)vs.'y 
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Fig. 4. 15, Br{B+ r;p+) vs. 7 
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Fig.4.16, Br{B° — > rip^) vs. 7 



FIG 4 



33 



) 50 100 150 200 250 300 350 

Fig.5.1, Air-{B+ r]'K+) 




50 100 150 200 250 300 350 

Fig.5.2, AfJ'iB" jj'K") 
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Fig.5.3, Aii,''{B+ riK'+) 
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Fig.5.4, Ai'^{B° TjK'°) 
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Fig.5.5, ^ V7r+) 
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Fig.5.6, Air'{B+ ri'K' + ) 
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Fig.5.7, A%''(B° ifK'°) 
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Fig.5.8, Ai'„''{B+ r,'p+) 
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Fig.5.9, Aii,^{B+ riK+) 
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Fig.5.10, A'^'(B° riK^') 
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